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The collision- free Boltzmann equation is used in the context of brane-J(-R) gravity to derive 
the virial theorem. It is shown that the virial mass is proportional to certain geometrical terms 
appearing in the Einstein field equations and contributes to gravitational energy and that such a 
geometric mass can be attributed to the virial mass discrepancy in a cluster of galaxies. In addition, 
the galaxy rotation curves are studied by utilizing the concept of conformal symmetry and notion 
of conformal Killing symmetry. The field equations may then be obtained in an exact parametric 
form in terms of the parameter representing the conformal factor. This provides the possibility of 
studying the behavior of the angular velocity of a test particle moving in a stable circular orbit. The 
tangential velocity can be derived as a function of the conformal factor and integration constants, 
resulting in a constant value at large radial distances. Relevant phenomenon such as the deflection 
of light passing through a region where the rotation curves are flat and the radar echo delay are also 
studied. 



PACS numbers: 04.020.-q, 04.50.Kd 



I. INTRODUCTION 

The question of dark matter is presently one of the 
most pressing open problems in cosmology. The galaxy 
rotation curves and mass discrepancy in a cluster of 
galaxies are two prominent observational evidence for 
the existence of dark matter. According to Newtonian 
gravity, galaxy rotation curves give the velocity of mat- 
ter rotating in a spiral disk as a function of the distance 
from the center of a galaxy; it increases linearly within 
the galaxy and drops off as the square root of 1/r out- 
side the galaxy. However, observations show that this 
is not the case and the velocity remains approximately 
constant. This provides for the possible existence of a 
new invisible matter distributed around galaxies which 
is known as dark matter [l], 0| • 

The mass discrepancy of clusters as another evidence 
for the existence of dark matter can be understood when 
estimating the total mass of a cluster in two different 
ways; summing the individual member masses within the 
cluster leads to a total mass which we shall call M. Al- 
ternatively, the virial theorem applied to a cluster would 
yield an estimate of the cluster mass which we call My. 
As it turns out, My is nearly 20 — 30 times greater than 
M and this difference is known as the virial mass dis- 
crepancy [l], Hj]. The best way to deal with the above 
discrepancy, it seems, is to postulate dark matter. There 
are several candidates for dark matter which can be cat- 
egorized as baryonic or non-baryonic, or, as hot or cold 
according to their velocity at the time when galaxies were 
just starting to form. 

In spite of many efforts to postulate various forms 
of matter as dark matter, there is as yet no non- 
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gravitational evidence for it. Moreover, accelerator and 
reactor experiments do not support the scenarios in 
which dark matter emerges. Therefore, one may conclude 
that Einsteins's gravity may break down at the scale of 
galaxies @. Thus, to deal with the question of dark mat- 
ter, modifications to Einstein field equations have 
become a flourishing method in recent years. One such 
modification is the brane-/ (i?) gravity. The idea that our 
4-dimensional universe might be a 3-brane embedded in 
a higher dimensional bulk has its roots in string theory 
[6|. One of the most successful of such higher dimen- 
sional models is the Randall- Sundrum (RS) scenario. In 
the RS scenario, our 4-dimensional universe is considered 
as a brane in a 5-dimensional bulk with an AdS geometry 
0, H[ ■ The RS model has had a great success in explain- 
ing the hierarchy problem. Another noticeable effort in 
this direction is the model proposed by Dvali-Gabadadzc- 
Porrati (DGP) @. In this model the bulk space is as- 
sumed to be flat but there is an additional 4D-induced 
gravity term appearing in the action. A self accelerating 
phase at late times is predicted for the universe in this 
model. However, the solutions corresponding to the self- 
accelerating branch suffer from ghost instabilities while 
the normal branch admits solutions which are ghost-free 
[ljjjill ■ The DGP model predicts modifications to grav- 
ity at large distances while the RS-type models modify 
gravity at small scales. Brane- world scenarios have paved 
the way for a new understanding of dark matter [l3j . In 
a similar vein, theories of gravity in which the Einstein- 
Hilbert (EH) action is replaced with a generic function of 
R, the Ricci scalar, have been quite useful in dealing with 
the question of dark matter [lj, EH ■ ^ would therefore 
be of interest to consider a brane-/ (R) gravity scenario 
incorporating both the above ideas to deal with the open 
questions in gravity. 

In addition, there are a number of other motivations 
for such studies. For example, since brane world sce- 
narios with a generic f(R) in a bulk with a single extra 
dimension can be described by a real scalar field, brane- 
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f(R) models have been studied within the framework of 
scalar-tensor type theories. Such theories have been used 
to investigate the stabilization of the distance between 
the branes in the context of the RSI model and the 
problem of the cosmological constant [17j. On the other 
hand, brane-/(i?) type theories have been used to answer 
questions such as dark energy and the present accelerat- 



ing phase of the universe 18[ . The present study could 



be considered as an attempt to narrow the gap in our 
understanding of one of the problems facing us in the 
ever increasing complexities of the inner workings of our 
universe. 

Generally speaking, one may start with a generic ac- 
tion involving f(R) without a cosmological constant in a 
RS scenario [19j in conjunction with the virial theorem in 
an attempt to account for dark matter. In doing so, the 
virial theorem together with observational data relating 
to the velocity of each member of a cluster, provides an 
estimate of the mean density and leads to a prediction 
of the total possible mass. The virial theorem has also 
been used in model f(R) theories with a cosmological 
20, 2l| , within the context of brane- world sce- 
using both the metric [l5[ and Palatini 
In this paper we generalize the virial 
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formalisms [1^ 
theorem in a brane- /(i?) scenario using the collisionless 
Boltzmann equation. In this process, extra terms will 
appear in the generalized virial theorem originating from 
the modified action in the bulk. The extra terms can 
be interpreted as a geometric mass and attributed to the 
mass discrepancy in clusters. 

According to recent observations, the tangential veloc- 
ity of matter moving around the center of a galaxy tends 
to a constant value as one moves away from the center of 
the galaxy. Such rotation curves have been studied in the 
context of brane-world scenarios by using the concept of 
conformal symmetry Q. We present similar conclusions 
in a brane-/(i?) model by using the idea of conformal 
Killing symmetry. In this regard, the spacetime is as- 
sumed to have, in addition to being static and spheri- 
cally symmetric, a conformal symmetry. If the vector 
field £ is the generator of such conformal symmetry, then 
the spacetime metric h is mapped conformally onto itself 
along the trajectories of £ 
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with C being the Lie derivative and ip is the conformal 
factor. There is a systematic method of searching for 
exact solutions, initiated by Herrera and co-workers [24l - 
[26j , where the field equations can be obtained in an exact 
parametric form, with the conformal factor taken as a 
parameter. This would provide the possibility of study- 
ing the behavior of the angular velocity of a test particle 
moving in a stable circular orbit. The tangential velocity 
can be derived as a function of the conformal factor and 
some integration constants. At large radial distances, we 
obtain a constant value for the tangential velocity by tak- 
ing suitable integration constants. The bendin g o f light 
is another issue worth studying in this region [5l,l27| . The 



computation of the deflection of photons passing through 
a region of flat rotation curves and their time delay are 
useful tools for testing the alternative theories of gravity, 
an example of which is presented in this work. 



II. GENERALIZED VIRIAL THEOREM IN 
BRANE-/(7?) GRAVITY 

In order to derive the virial theorem for galaxy clusters, 
it is necessary to introduce the brane-/ (i?) model we are 
interested in. We also need to know the Boltzmann equa- 
tion governing the evolution of the distribution function 
in cluster of galaxies. By taking the cluster of galaxies as 
a system of identical and collisionless point particles, we 
utilize the relativistic Boltzmann equation together with 
the field equations to find the generalized virial theorem. 



A. The brane-/(i?) gravity 

In a brane-/(i?) model, the 5-dimensional bulk action 
is taken as 



S = J d 5 x^g~[f(R)+£ m ], 



(2) 



where C m is the matter lagrangian, g is the bulk metric 
and R is the bulk Ricci scalar [19| ■ Variation of S with 
respect to the bulk metric qab yields 



F{R)R AB - ^ 9 ABf{R) + 9abOF{R) 
- V A V B F(i?) =k 2 5 T ab , 

_ dftm 



(3) 



where F(R) — ai }^' ■ The effective Einstein field equa- 
tions in the bulk can be written as 



Gai 



where 



rptot 

1 AB — 



F(R) 



Rab - -^Rqab — Tab, 



k\T ab - (-RF(R) - -f(R) 



(4) 



aF(R))g AB + V A V B F(R) 



(5) 



Using the SMS procedure [28[ , the field equations on the 
brane are given by 

G^ v = SttGnt^ + k^tt^ + Q^ u - E^ u . (6) 

We note that 



Gn — 



48tt' 



(7) 



where A is the brane tension, t^ v is the energy momentum 
tensor on the brane and ir^ is defined in terms of as 



1 

12' 



a/3 



h.njT 1 

24 M 



(8) 
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The electric part of the Weyl tensor is given by 



E„u = C* CD n A n c h B h°, 



(9) 



where n is the unit vector normal to the 4-dimensional 
brane and Iiab — 9ab — riAriB is the induced metric on 
the brane. Furthermore, we have 



Q uu 



n(i?)v 

+ n A n B h^) 



2y a VbF(R). uAuB 



F(R) 



(10) 



y=o 



and 



n(i?) 



4 DF(R) 
15 F{R) 



10 V 2 



F(R)j 



lf(R) \uF{R), 



(11) 



where y is the extra dimension and the brane is located 
at y = 0. Now, if we define <2 M „ — P M „ as 



we obtain the field equations 



G 



//// 



(12) 



(13) 



Using t-; = diag(-p e ff,pr ff ,P^- ff ,P^ ff ) and ti"t 
— 1, the field equations become 



-x 



r 

jr \2 



-^)--^2= -^GNPeff + ^ 



X \2 



+ (p; // )^-2(p; f/ )(p^ / ) + (p^ / ) 



eff' 



(Peff) 2 

8itGnPx, 

(18) 



+ 2(P e L // )( Pe// )-(PJ // ) 2 + (P e L // ) 2 



(Pe//) 2 

SttGatPJ,, 



(19) 



and 



v' X' v'\' v" v n , 



^[(Pe//) 2 + (p e //)(PJ // ) + (P e 7 / ) 2 



+ (Pe//)(P e L // ) - (P e 7 / )(P e L // )] + 8nG N P^. (20) 

We note that the (90) and (</9ty?) components of the field 
equations are similar. 



where J^v may act as a new matter source on the brane 
induced by the f(R) action in the bulk. It is convenient 
to represent this new matter by 



7? = (-px,Pi,P£,P£). 



(14) 



Let us now consider an isolated and spherically symmet- 
ric cluster described by the metric 



ds 2 



^ {r) dt 2 + e x{ - r) dr 2 + r 2 d0 2 + r 2 sin 2 Odip 2 , (15) 



living on the brane. Suppose that the clusters are con- 
structed from galaxies which are acting as identical and 
collisionless particles and described by the distribution 
function fs ■ The energy momentum tensor may be writ- 
ten in terms of Jb as [29[ 



fBmUftiivdu, 



(16) 



where m is the cluster's member mass, u is the four ve- 
locity of the galaxy and du = dUrd ^ du ^ [ s th e invariant 
volume element of the velocity space. The energy mo- 
mentum tensor of the matter in a cluster is given by an 
effective density p e // and an effective anisotropic pres- 
sure, with radial p r e ff and tangential P^ff components 
201. In other words, we have 



Peff 



P{u 2 t ) 



P 



eff 



p(u 2 }, 



P e A ff =p(u 2 e )=p(ul). 



(17) 



B. The generalized virial theorem 

To derive the virial theorem, we need the Boltzmann 
equation which governs the evolution of the distribution 
function. By integrating this equation on the velocity 
space and using the gravitational field equations, the 
virial theorem can be obtained. We consider the cluster 
as an isolated spherically symmetric system described by 
equation (fT"5|). Furthermore, we assume that the galaxies 
in the cluster behave like identical, collisionless point par- 
ticles. The distribution function is denoted by fg which 
obeys the general relativistic Boltzmann equation. 

In many applications, it is convenient to work in an 
appropriate orthonormal frame or tetrad eJJ(x), where 
<7 My e°e^, = rj ah . In the case of the spherically symmetric 
line element given by equation (1151), w e introduce the 
frame of orthonormal vectors as [15|, [29[ 



o u/2 5 o 



rSf. , ef = r sin 65,, . 



(21) 



We also note that the tetrad components of u M can be 
written as u a = it M e^. In tetrad components, the rela- 
tivistic Boltzmann equation can be written as 



afl dfs 



i b c ®fs 



0. 



(22) 



where Jb — /b^^ 1 , - "") is the distribution function and 
Ibc = e 1j,-v e b e c are t ne Ricci rotation coefficients [H, [20|, 
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|29[ . We may assume that Jb depends only on the radial 
coordinate r. Thus, the relativistic Boltzmann equation 
becomes 



dfB_ 
dr 

dfa 
du 3 



1 2 dv II 



-u 3 



2 -r uj \ df B _ l u ( df B 
r J dui r \ du 2 

^u 3 e x / 2 cote(u 2 ^--u 3 ^ ]=(). 
r V ^"3 <ju 2 J 



df B \ 

i I 



(23) 



The term proportional to cot 9 must be zero since we have 
assumed the system to be spherically symmetric. Let us 
take 



acting as a new matter source on the brane. It carries 
all the effects induced on the brane by the f(R) action 
in the bulk. We therefore have 



2v' v'A' „ v . 
— + v " + — = 8nG N p + 8irG N p x . 



—X 



(29) 

To move on, we assume that v and A' are slowly vary- 
ing, i.e. v 1 and A' are small and their product can be 
neglected, and e~ A « 1 inside the cluster. We then have 



u = u t . m = u r , u 2 = u e , u 3 =u 



(24) 



Multiplying equation (J23J) by mu r du, integrating over the 
velocity space and assuming that fg vanishes sufficiently 
rapidly as the velocities tend to ±oo, we find 



WD] 



1 dv 
2dr~ 



p[{i 
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Now, it is helpful to multiply equation (|25|) by Airr 2 and 
integrate over the cluster volume to obtain 



l - P [{ul) + «)] + - r p{ul) = 0. (25) 



u r .) + (u, 

R 



\j\ 4nr 2 dr 



-\[ P [<«?> + %^ dr = °> (26) 

where R is the radius of the cluster. A useful equation 
is obtained by summing all none zero components of the 
field equations (I18H20|) 



— A 



2v' v'X' 



— --i- + v" + -T)= ^G N p[(u 2 ) + (u 2 ) 



+ («?) + «)] +^f[2<^) 2 + 2<«?><«|> 
+ (u 2 t ){u 2 )] + 87rG N {p x +P r x + 2Pi], 



(27) 



where we have used equation (|17[) . Using 



u 



(u. 



(u£) + (ibg) + (u^) and assuming that the galaxies in the 
cluster have velocities much smaller than the velocity of 

A) 



light [I5j,l22j, that is (u 2 r ) a (u 2 ) « (u 2 ) < (u») » 1, we 



obtain 



= SttGnP 



2v' v'W „ i 

— + v" + - 

r 2 

I/^+SttGjv + 2P^) 



(28) 



Since for clusters of galaxies the ratio of the matter den- 
sity to the brane tension is much smaller than unity, 
j- <C 1, we can neglect the quadratic term in the mat- 
ter density in equation (|2"5|). We may also define px = 
{fix + Px + 2-Pjf), where is a pure geometric term 



On the other hand, using the above assumptions, one 
may write equation (|2"6")l as 



1 



dv 



2K- - Airr 6 p—dr = 0, 
2 Jo or 



(31) 



where 



K 



2tt P [(u 2 r ) + (u 2 ) + (ul)] r 2 dr, (32) 



is the total kinetic energy of the galaxies. Multiplying 
equation (|30[) by r 2 and integrating yields 

1 dv 

G N M(r) = -r 2 ^-- G N M x (r), (33) 

where we have used M = /„ dM(r) = / J Anpr 2 dr 
as the baryonic mass. We have also defined Mx = 
4irpxr 2 dr as the geometric mass of the system. Now 
consider the following definitions 



G N M(r) 



dM(r), 



(34) 



ax= f R ^M x (r) dM{ry 



and 



Multiplying equation (133)) by dM ( r ) which is equal to 
47rpr dr and integrating gives 



1 



dv 



= ^ - o / ^p—dr 



(36) 



where is the gravitational potential energy of the sys- 
tem. Finally, using equation (|31l) leads to the generalized 
virial theorem 



2K + n - n x = o 



(37) 



Alternatively, the above equation can be written in the 
form 



2K-G 



N 



Jo r Jo 



dM 
r 



(38) 
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For convenience we introduce the radii Ry and Rx 



and 



Rv = 



R x — 



M 2 



fR Mir) ,„ w \ 

Jo -^dMir) 



Ml 



Furthermore, the virial mass My is defined as [if 



2K 



Rv 



(39) 



(40) 



(41) 



Substitution of these definitions in equation (I3"51) leads to 



My 

M 



M\R V 
M 2 R 



x 



(42) 



From observation, the relation Mv /M > 3 holds true for 
most of the galactic clusters. Therefore one can easily 
approximate the last equation 



My 

M 



M\Ry 

M 2 R X ' 



(43) 



Now one may note that some geometric terms appear- 
ing in the Einstein field equations above could effectively 
play a role in the gravitational energy. These geomet- 
ric terms may be attributed to a geometric mass at the 
galactic or extra galactic levels and may be interpreted 
as dark matter. On the other hand, dark matter is the 
main contribution of mass in clusters. It means that the 
contribution of the baryonic mass is negligible in compar- 
ison with dark matter. The total mass of the cluster can 
then be estimated as M to t ~ Mx • We also know that the 
virial mass is mainly determined by the geometric mass. 
It means that the geometric mass could be a potential 
candidate for the virial mass discrepancy in clusters. As 
a result, we conclude that 



M x « My a M to f 
Therefore, equation (|43p can be written as 

-Rx 



My w M 



Ry 



(44) 



(45) 



This shows that the virial mass is proportional to the 
normal baryonic mass in a cluster and the proportionality 
constant has geometrical origins. 



III. VACUUM SOLUTION 

In this section, we consider the vacuum solution of the 
theory. Assuming r M „ = 0, the field equations (jl8H20[) 
reduce to 



A' 



1 



= %-kGnpx, 



(46) 



8ttG n P x , 



(47) 



and 



v 
2~7 



27 



AV 



v 

T 



8irG N Px-. (48) 



Since the above system of equations contain five unknown 
quantities, it is under determined. To obtain the un- 
known quantities we require two extra relations. There- 
fore, it is convenient to consider 



Px = P x 



X > 



(49) 



as an equation of state. We also assume that the space- 
time admits a conformal group of symmetries Q. In 
other words, if the vector field £ is the generator of such 
a conformal symmetry, then the spacetime metric h is 
mapped conformally onto itself along the trajectories of 
t 



C e h 



(50) 



where C is the Lie derivative operator and ip is the con- 
formal factor [24Tj26| . It should be noted that such extra 
relations affect the form of f(R) through relation (|49l) . 
Let us consider a general form for the vector field £ 



d_ 

dip 
(51) 



Substituting this conformal vector in equation (|50|) and 
assuming tp = ip( r )i one obtains 

(52) 

where F(ip) and G{9) are arbitrary functions and k is an 
integration constant [30( ] - This procedure also leads to 
the metric components in terms of the conformal factor 



Mr) = c 2 r 2 exp f -2k I — 

rip 



and 



X[r) _ 



(53) 



(54) 



where B and C are integration constants [3(|. Substitu- 
tion of equations (l53l) and (154)) into field equations ((46])- 
leads to 



2i>' 
rip 



1 



1 



lP 2 

B 2 



2k 

r 2 ip 



2 = 8itGnP 



N^X, 



(55) 



(56) 
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and 



f_ / 1_ k 2 

£2 ' ° 



2k 2i>' 
r 2 ip 2 r 2 ip rip 



= 8ttG n P^. 



(57) 



Using the equation of state (|49[) and equating equations 
(155]) and ([5d| we obtain 



r^' + 2?A 2 - kip - B 2 = 0. (58) 
This equation can be solved to give r in terms of ip 



r 4 = 7? 4 - 

°|2^ 2 -/cV-S 2 |' 

where i?o is an integration constant and 



G(ip) = exp 



-2fc 



VA: 2 + 8B 2 



: tanh 1 



-Aip + fe 



Vfc 2 + 8B 2 



(59) 



(60) 



Now, by substituting all quantities in terms of ip, the 
energy density and pressure of what may now be called 
the X-matter can be written in terms of the conformal 
factor as 



8TrG N px = 8irG N P r x 



^-2ip 2 + kiP + B 2 
B 2 R 2 ^Gm~ 



(3V> 2 - 2k4> - B 2 ) 



(61) 



and 



^J-2iP 2 + kijj + B 2 

b 2 r 2 ^gWT~ 



(3^ 2 



2B 2 - k 2 



(62) 



The sum of equations ffBTT) and (p2]l is another useful 
relation 



8irG N (p x +P^) = 



1 



B 2 r 2 



{k 2 -2ktjj + B 2 ), (63) 



which can be used to define the equation of state at large 
distances. It is clear that in the limit r — > oo, the X- 
matter has an equation of state px = Px = ~^x- 

It should be emphasized that to find the vacuum solu- 
tion we only ignored ordinary matter on the brane. The 
presence of f(R) in the bulk action results in modifica- 
tions which can be traced to /(i?)-nonlinearity effects 
emanating from and E^. Such effects constitute 
what is known as the X-matter. The vacuum solutions 
obtained here may thus be usefully utilized to account for 
the flat rotation curves in a region where the contribution 
of X-matter is dominant. 



IV. GALAXY ROTATION CURVES IN A 
CONFORMALLY SYMMETRIC SPACETIME 

The observational data show that the rotational veloc- 
ity increases linearly within the galaxy and approaches 



a constant value of about 200km/ s as one moves away 
from the center [2;j. In this section, we shall consider a 
test particle which moves in a circular time-like geodesic 
orbit and study its tangential velocity. The Lagrangian 
of the system is given by 



2L 



"W 



dtY 



Mr) 



Tr ' ' (64) 



where r is the affine parameter along the geodesic and 
dQ 2 = d6 2 + sin 2 9d(p 2 . From the above Lagrangian one 
obtains 



E = e' y t = const. 



and 



r 2 sin 2 i 



const., 
(65) 



where E and 1$ are the energy and 0-component of the 
angular momentum of the test particle respectively which 
are conserved quantities and an over dot represents differ- 
entiation with respect to t. Although the ^-component 
of the angular momentum is not a constant of motion, 
the total angular momentum, I 2 = I 2 , + (l^j sin#) 2 , is a 
conserved quantity which can be written as I 2 — r 4 f2 2 
}3~fl ] . Equation of the geodesic orbit is given by 



r 2 + V(r) = 0, 



where V(r) is the potential 
V(r)-- 



-A 



I 2 



(66) 



(67) 



Let us now study the motion associated with stable cir- 
cular orbits. This kind of motion can be obtained by the 
conditions 







and 



dV 
dr 



0. 



(68) 



The potential obtained subject to the above conditions 
describes an extremum of the motion. To have a mini- 
mum, the condition 



d 2 V 
dr 2 



>0, 



(69) 



is also required. Using these conditions, one obtains the 
energy and total angular momentum [3l|, H1J 



E 2 



2e v 



and 



I 2 



(70) 



For an inertial observer far from the source who measures 
the spatial components of the velocity, normalized to the 
speed of light, the line element can be rewritten as ds 2 = 
-dt 2 {l-v 2 ) ;31;j, where 



v 2 = e 



+ r* 



dnv 
~dj) 



(71) 
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Using the condition for circular orbits, the tangential ve- 
locity can be expressed as 



u tg 



r 2 e 



dttY 
~dt) 



(72) 



The tangential velocity can be rewritten in terms of the 
conserved quantities 



s" I 2 



E 2 ' 



or in the following alternative form 



u tg 



TV 
~2~ 



(73) 



(74) 



where equation ((70)) is used. This relation shows explic- 
itly that the tangential velocity depends on the time-time 
component of the metric only. Previously, we obtained 
the metric components in terms of the conformal factor 
in equations ([53")) and ([541 . Therefore, using the relation 
between v and ip, the angular velocity in terms of ip is 
given by 



(75) 



Using equation ([75j) and (|54|) . the metric coefficient 
exp(A) can also be expressed in terms of the angular ve- 
locity 



exp(A) 



B 2 

k 2 



(1 



'I? 



(76) 



Although, we cannot write vt g as a function of r explic- 
itly, it is possible to find Vt g in some regions. For example, 
we can derive the value of tangential velocity at infinity. 
It is clear from equation (1591) that we have r — > oo for 
"0 — "01,2, where 



1>i 



± 



k 2 
16 



IP 
~2~' 



(77) 



Therefore one may obtain the tangential velocity at in- 
finity from equation (|75|) as follows 



Vtgc 



\ 



5 ± 



ML _i_ Ml 

16 ' 2 



(78) 



Knowing the value of the tangential velocity at infinity 
from observations one can derive the relation between k 
and B. The variation of v tg as a function of r, using equa- 
tions (|75l l58 j) . for vtgoo = 200km/ s is shown in FigfTJ It 
should be noted that parameter Rq is an integration con- 
stant from equation (I5"9")) and helps to re-scale the curve 
only. Furthermore, B should not be looked upon as a 
fine tuning parameter as it can take values not necessar- 
ily close to 1. The particular choices for B in FigfT]are 
only made so that the curves are better distinguished. 



v tg (km/s) 
200 : 

190 ; 

180 r 

170 ; 
160 ; 
150 ; 

140 - 




FIG. 1: Variation of the tangential velocity as a function of 
r/Ro, where the dotted line represents B = 1.1000007, the 
dashed line B = 1.1 and the solid line B = 1.0999998. 



Now, it is convenient to obtain geometric mass of the 
galaxy and its dependence on r. If we integrate equation 
(1461). we will find 



= 1 - 



BttG 



pxr 2 dr. 



(79) 



Using definition Mx — 47r px{r')r l2 dr' , we also obtain 

2Gjv Mx 



= 1 - 



(80) 



One may now introduce the radius r — R for the vacuum 
boundary where the contribution of the baryonic mass 
vanishes, pb ~ 0. Of course, the continuity of the met- 
ric coefficients e A and e v across the vacuum boundary 
of the galaxy should be preserved. For simplicity we as- 
sume that inside the baryonic matter with density pb the 
nonlocal effects of the bulk can be neglected. Therefore 
at the vacuum boundary the metric coefficients can be 
written as 



_ A 2G N M b 
R ' 



(81) 



where M& = 47T J q pb{r')r' 2 dr' is the total baryonic mass 



of the galaxy. By equating equations (ffrJ]) and (EI 
r = R, one has 



at 



*1 
B 2 



1 



2G N M, 
R 



1 



J tg) 



(82) 



Using the above formula, it is possible to find a relation 
between Mx and Mb from equations (|76|) and ([80]) as 
follows 



M X = ^r, 
R ' 



(83) 



showing that Mx varies linearly with r, consistent with 
observations. 
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V. THE LAGRANGIAN 

To find the form of the Lagrangian f(R) of the theory 
on the brane, we need the metric. Here, we concentrate 
on the flat rotation curve region for which v tg is constant. 
According to equation ([75)1. we have V J = k/(l — uf„). 
Therefore, the constancy of v tg in the region of flat rota- 
tion curves leads to the constancy of ip. As a result we 
may write 



C 2 r 2 < 



(84) 



and 



So far we have found the metric in the region of flat 
rotation curve as given by In this region we use the 
field equations © on the brane to determine the form 
of the Lagrangian. By neglecting the contribution of the 
baryonic matter, the field equations reduce to 



IV Fu,Ui 



where and are given by equations (0) and (fTU)) 
respectively. To solve the field equations, we start with 
the bulk metric and write it as 



(85) 



ds 1 



-N{y)r a dt 2 + N{y)Adr 2 + r 2 dil 2 + dy 2 , (89) 



obtained from equations ([5U]) and (|54"1) . Therefore the 
form of the line element can be written as 

2 B 2 

ds 2 = -C 2 r 2v ^dt 2 + — (1 - v 2 ig ) 2 dr 2 + rW, (86) 

K 

where the relation between k and B is known according 
to equation ([52)1 and C can be defined according to the 
continuity of the metric coefficient e v across the vacuum 
boundary of the galaxy. Using equations (|5Tj) and (|8"4"|) 
we have 



C 2 = R- 2v *s 1 - 



2G N M b 
R 



(87) 



where A is a constant and N(y) is a function of the extra 
dimension. To obtain a consistent line element on the 



brane as in 
and N(Q) - 
becomes 



we assume a — 2v\ L, A — 



C 2 



Cik? v tg) 2 

With these assumptions, the metric 



ds 2 



-C 2 r a dt 2 



Pdr 2 



l dVL 2 



(90) 



where /3 = C 2 A. To continue, we can find f{R) from 
the trace of equation ([J) and substitute it in the brane 
field equations. The non zero components of the field 
equations become 



2 , ,„ . 2 , 1 IN 2 7N \ 2 F r f4 \ 4 F rr 1 ( „z . NF 



(-4/3 + 4 -a z -a) + —, _-_ +_-£ a + _£1 2F + — , = 0, (91) 



15/3r 2 ' r ' 5C 2 C 2 3 5/3r F \3 / 15/3 F 5F I 3C 2 



:( _ w+4 + 9lI _^ ) + i £ ™ +7 4%«» + 2,-A% i 2 ^ + ££ =0 , (92) 



15^r 2 v ^ ; 5C* 2 \ C 2 3 I 15/3r F 5/3 F 5F I 3C 2 



* rtf + to »_ _ 1) 1 U^p-D + ^-i -0, (93, 



15/3r 2V ; 5C 2 I 4C 2 3 / 15/3r F ; 15/3 F 5F I 3C 2 

where F = F(r,y), F r = §E | y=0 and F rr = |„ =0 . In addition, we have defined F = f| | y=0 and F = ^=0, 
with the same notation for N = N(y). We emphasize that N and its derivatives are constant on the brane. 



From equations (1911) and (|92p we find where C\ and C2 are arbitrary integration constants and 

ni,2 = (a + 1 ± y/a 2 + 10a + l)/2. We limit ourselves 
—r 2 F rr + arF r + 2aF = 0, (94) to the monotonically decreasing solution so that F(r) = 

C^r" where n — n 2 < [Hj]. Now, we can derive F and 

with general solution 

F(r) = Cir" 1 +C 2 r' 12 , (95) 
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F from equations (|9lj) and (|93[) with the result 

7i , 37V /v\ 



At* 



F(r) - C 2 t 



r 2 1 4C 2 TV / ' 



(96) 



where /i = §^hVa 2 + 10a + 1(1 - a) + 3 - 3a 2 - 2/3] 
and 



(97) 



with / 2 = ■X [- Va 2 + 10a + 1 + 3a - 2/3 + 3] . Using the 
above formulae and the trace of equation Q one finds 



f(R(r,y = 0)) = C 2 r r ' 



h 



2N 2 



(98) 



where / 3 = fhVa 2 + 10a + 1 - /3 - a 2 + 2a + 2]. 

At y = 0, one may use the bulk metric to find the Ricci 
scalar as 




2C 2 



2N 



1). (99) 



It is clear that one may substitute r in terms of R in 
equation (198)) to derive the action on the brane. How- 
ever, the same is not possible for the bulk action. This 
is due to the fact that we have derived all equations on 
the brane by setting y = and therefore there remains 
no information on the extra dimensional dependency. As 
a result, starting from a given metric and using the ob- 
served physical properties of a particular galaxy such as 
Vt g , we may derive the form of f(R) on the brane for the 
region of flat rotation curves of the given galaxy. Such 
a f(R) function may now be used for the fiat rotation 
curves of any other galaxy to find the physical parame- 
ters and hence the metric. 



VI. LIGHT DEFLECTION ANGLE 

To calculate the light deflection angle we use the metric 
found in the region of flat rotation curves, in the form 
given by (156")) . The bending of light results in a deflection 
angle A<fi given by 



A</> = 2 | <t>{r c ) - 0oo | 



(100) 



where 4>oo is the incident direction and r c is the coordi- 
nate radius of the closest approach to the center of the 
galaxy [39]. Using the geodesic equations, deflection of 
light can be written as 



(r c ) - 0c 



Y(r) 



v{r a )-v{r) 



dr 
r 

(101) 



8.xl(T 6 



6.x 10 



4. x 1(T 6 



2.xl(T 6 




r c /r<, 



FIG. 2: Variation of the light deflection angle as a function 
of r c /ro for v tg = 300km/ s. 



We consider the case where r c is in the region of flat 
rotation curves. Therefor, the integral is split into two 
parts 



r c ) - 0oo = 

ro n 



r x Of. -2) 



-1/2 



dr 
r 



2G N M 



1 - 



2G N M / x 2 
r 



2G N M 



-1/2 



dr 
r 



(102) 



where ro is the radius for which the contribution of 
dark matter vanishes. The first integral is attributed to 
the region of flat rotation curves with the metric given 
by equation (1861) where dark matter is dominant. The 
second integral relates to the exterior region with the 
Schwarzschild metric ^3. As a result we obtain 



A0 : 



i-^K 2 9 -i) 



2 arcsin 



'•() 



arctan 



2G N M 



1 



(S) 



OL-2) 




ro 



ro 



(103) 



as the total light deflection angle. It is clear that in the 
case R = rg — r c , the deflection angle only originates 
from baryonic matter. In this particular case, we have 
A0 = AGnMi,/R — Ax 10~ 6 which is consistent with the 
general relativity prediction. The effects of dark matter 
on the deflection angle is shown in Figf2] where variation 
of At/> as a function of r c /ro is plotted. It is necessary to 
point out that this result is valid for the region R < r c < 
ro. In general, the point of the closest approach, r c , can 
be anywhere. 
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VII. RADAR ECHO DELAY 

The time needed for light to travel by a massive celes- 
tial body is longer than the same time calculated by New- 
tonian gravity. This effect originates from the existence 
of gravitational field in the vicinity of these massive ob- 
jects and can be detected using a radar echo delay. Such 
a delay was measured by a method suggested by Shapiro 
[ioj . The time needed for a photon to propagate from 
one point with r = r\, (f> = <f>i, 6 = ^ to another point 
with r = r2, 4> = 4>2, = ^ can be obtained by use of 
geodesic equations 39]. The elapsed time for a photon 
traveling between r c and r is given by 



t(r, r c ) 



r e A(r)/2 



"(rc) \ r ) 



dr, (104) 



where r c is the closest approach to the celestial body. 

To find the time delay we suppose that r is outside 
the region formed by the galaxy and r c is in the region 
of flat rotation curves. Thus, we split the integral (|104|) 
into two parts according to 



t{r,r c ) = 
rr B 
C~k 



1 - 



1 

2G N M 



r 2 < 1 



1 - 



r 
r c 

2G N M 



(2< 9 -2) N 



-1/2 



dr 



Y 2G jy M y j. 



-1/2 



dr 
r 



(105) 



where we can integrate the first term exactly and the 
second term by using the Robertson expansion [39| to 
find 



t(r, r c ) 



B 



1- I 1° 



Ckr 



2uf„-2 



2G N M In ( r + ) + G N M ( "—^ 




1/2 



(106) 



The term yr 



is what we would expect if light trav- 



eled in a straight line with unit velocity. The first term 
in the above expression is the effect of the flat rotation 
curves region and the other terms represent the effect of 
general relativity. It is clear that in the absence of the 



region containing the flat rotation curves, r c = tq, one 
gets the time delay for general relativity only. 



VIII. CONCLUSIONS 

In this work we have considered a brane-/(i?) model as 
a possible candidate to explain the question of dark mat- 
ter. This was achieved by considering the field equations 
on the brane obtained through the projection of the field 
equations in the bulk together with the assumption that 
our spacetime admits a family of conformal symmetries. 
The additional terms thus appearing in this process in 
the field equations can be considered as a source for dark 
matter. We used this procedure to explain the problem 
of what is known as the flat rotation curves and obtained 
an expression for the tangential velocity of a test particle 
moving in such a region, consistent with present obser- 
vational data. In addition, the angle representing the 
deflection of light and the time representing the radar 
echo delay of photons passing through such a region were 
calculated. 

To look further afield, one may be interested to in- 
vestigate the stability of the model and its consistency 
with Parameterized Post-Newtonian (PPN) parameters. 
Four dimensional f(R) gravity models have been con- 
sidered before in order to study their stability [H, H3| 
and PPN parameters [H, |36j]. This is basically done 
by investigating weak field limit of the theory involved. 
The same route can be taken here, leading to the pos- 
sibility of direct comparison of predictions of the model 
with observational data. The weak field limit would also 
afford the possibility of obtaining the potential energy 
which is needed to find the Tully-Fisher relation [33, [38| ■ 
This important relation establishes a connection between 
the rotational velocity of a spiral galaxy and its luminos- 
ity. Since the luminosity is proportional to the mass of 
the galaxy, the Tully-Fisher relation connects the rota- 
tional velocity of a galaxy with its mass. In addition, 
the virial theorem provides a relation between the ki- 
netic energy (velocity) and potential energy (relating to 
mass or luminosity) in self-gravitating systems in equi- 
librium |4lj . As a result, the potential energy found from 
the weak field limit equations accompanied by the New- 
tonian virial theorem [37| points to a possible derivation 
of the Tully-Fisher relation in the present context which 
furnishes another possibility for investigating the validity 
of the model. 



[1] M. Persic, P. Salucci and F. Stel, Month. Not. R. Astron. 
Soc. 281 (1996) 27, |arXiv:astro-ph/ 9506004 
A. Borriello and P. Salucci, Month. Not. R. Astron. Soc. 
323 (2001) 285, |arXiv:astro-ph/000Tu82| 
P. Salucci, A. Lapi, C. Tonini, G. Gentile, I. Yegorova 



an d U. Klein, Mon. N ot. Roy. Astron. Soc. 378 (2007) 

41, |arXiv:astro-ph/0703115| 
[2] J. Binny and S. Tremaine, Galactic Dynamics, Princeton 

University Press, Princeton, 1987. 
[3] M. K. Mak and T. Harko, Phys. rev. D 70 (2004) 024010, 



11 



|arXiv:gr-qc70 404104 
[4] M. Milgrom, Astrophys. J. 270 (1983) 365; 
R. H. Sanders, Astrophys. J. 136 (1984) L21; 
J. Bekenstein and M. Milgrom, Astrophys. J. 286 (1984) 
7; 

R. H. Sanders, Astron. Astrophys. 154 (1986) 135; 

P. D. Mannheim, Astrophys. J. 419 (1993) 150, 

|arXiv:gr-qc/9407010| 

J. W. Moffat and I. Y. Sokol ov, Phys. Lett. B 378 (1996) 
59, |arXiv:astro-ph/9509143| 

P. D. Mannheim, ^Astrophys. J. 479 (1997) 659, 



arXiv:astro-ph/9605085 
M Milgrom, 



Rev. 46 (2002) 741, 



New Astron. 
|arXiv:astro-ph/0207231| 
M. Milgrom and R. H. Sander s, Astrophys. J. 599 (2003) 
L25, [arXiv:astro-ph/0309617| 

J. D. Bekenstein, Phy s. Rev. D 70 (2004) 083509, 
|arXiv:astro-ph/0403694| 

M. D. Roberts, Gen. Rel. Grav. 36 (2004) 2423; 

J. R. Brow nstein and J. W. M offat, Astrophys. J. 636 

(2006) 721, |arXiv:astro-ph/0506370| 

J. R. Brownstein and J. W. Moffat, Mo n. Not. Roy . As- 
tron. Soc. 367 (2006) 527, |arXiv:astro-ph/0507222| 
C. G. Bohmer and T. Harko, Class. Quantum Grav. 24 

(2007) 3191. larXiv:0705. 24961 

[5] T. Harko and K. S. Che ng, Astrophys. J. 636 (2005) 8, 
|arXiv:astro-ph/0509576l 

[6] M. Green, J. Schwarz and E. Witten, Super-string The- 
ory, Cambridge Monographs on Math- ematical Physics, 
two volumes, Cambridge University Press, 1987; 
Lust and Theisen, Lectures on String Theory, Springer 
Verlag, 1989; 

J. Polchinski, String Theory, Cambridge Monographs on 
Mathematical Physics, two volumes, Cambridge Univer- 
sity Press, 2004; 

B. Zweibach, A First Course in String Theory, Cam- 
bridge University Press, 2004. 
[7] L. Randall a nd R. Sundrum Phys. Rev. Lett. 83 (1999) 

3370, |arXiv:hep-ph/9905221| 
[8] L. R andall and R. Sundrum Phys. Rev. Lett. 83 (1999) 

4690, |arXiv:hep-th /9906064 
[9] G. Dvali, G. G abadadze and M. Porra ti, Phys. Lett. B 

485 (2000) 208, |arXiv:hep-th/00050"l6| 
[10] A. Nicolis and R. Ra ttazzi, JHEP 0406 (2004) 059, 

|arXiv:hep-th/0404159| 
[11] D. Gorbunov, K. Koyama and S. Sibiryakov, Phys. Rev. 

D 73 (2006) 044016, |arXiv:hep-th/0512097| 
[12] C. Charmousis, R. Gregory, N. Kaloper and A. Padilla, 

JHEP 0610 (2006) 066, |arXiv:hep-th/06040"86| 
[13] N. Dadhich, R. Maartens, P. Papadopoulos and V. Reza- 

nia, Phys. Lett. B 487 (2000) 1, |arXiv:hep-th/0003061~j 
[14] A. S. Sefiedgar, K. Atazadeh and H. R. Sepangi, Phys. 

Rev. D 80 (2009) 064010, 1 arXiv: 0908. 26931 
[15] C. G. Bohm er, T. Harko and F. S. N. Lobo, JCAP 0803 

(2008) 024 , arXiv:0710.0966l 

[16] W. D. Gold berger and M. B. Wise, Phys. Rev. Lett. 83 
(1999) 4922, |arXiv:hep-ph/9907447l 

O. DeWolfe, D. Z. Freedman, S. S. Gubser and A. Karch, 
Phys. Rev. D 62 (2000) 046008, arXiv:hepth/9909134. 
[17] N. Arkani-Hamed, S. Dimopoulos, N. Kaloper 
and R. Sundrum, P hys. Lett. B 480 (2000) 193, 
|arXiv:hep-th/0001197l 

S. Kachru, M. Schulz and E. Silverstein, Phys. Rev. D 
62 (2000) 045021, [arXiv:hep-th/0001206| 



[18] 

[19] 

[20] 
[21] 



K. Atazadeh, M. farhoudi and H. R. Sepangi, Phys. Lett. 
B 660 (2008) 275, larXiv:080Tl 398 

A. Borzou, H. R. Sepangi, S. Shahidi and R. Yousefi, 
EPL, 88 (2009) 29001. larXiv:0910.1933l 
J. C. Jackson, Month. Not. R. Astr. Soc. 148 (1970) 249. 
M. Nowakowski, J. C. Sanabria and A. Garcia, Phys. 
Rev. D 66 (2002) 023003, |arXiv:astro-ph/0105212] 
A. Balaguera-Antoh'nez, (J. GT Bohmer and M. 
Nowakowski, Class. Quantum Grav. 23 (2006) 485, 



arXiv:gr-qc/0511057 



[22] 
[23] 
[24] 
[25] 
[26] 
[27] 
[28] 
[29] 



A. Balaguera-Antoh'nez, D. F. Mota and M. Nowakowski, 
Mon. Not. Ro y. Astron. Soc. 382 (2007) 621, 
larXiv:0708T2980l 

T. Har ko and K. S. Ch eng, Phys. Rev. D 76 (2007) 
044013, larXiv:0707.lT28l 

Malihe Heydari-Fard and Mohaddese Heydari-Fard, 
Phys. Rev. D 84 (2011) 024040. 

L. Herrera, J. Jimenez, L. Leal, J. Ponce de Leon, M. 
Esculpi and V. Galina, J. Math. Phys. 25 (1984) 3274. 
L. Herrera and J. Ponce de Leon, J. Math. Phys. 26 
(1985) 2302. 

R. Maartens, D. P. Mason and M. Tsamparlis, J. Math. 
Phys. 27 (1986) 2987. 

U. Nucamendi, M. Salgado and D. Sudarsky, Phys. Rev. 

D 63 (2001) 125016, |arXiv:gr-qc/00lT049"| 

T. Shiromizu, K. Maeda, M. Sasak iTThys. Rev. D 62 

(2000) 024012, |arXiv:gr-qc/9910076| 

R. W. Lindquist, Annals of Physics 37 (1966) 487; 

R. Maartens and S. D. Maharaj, J. Math. Phys. 26 (1985) 

2869; 

S. Bildhauer, Class. Quantum Grav. 6 (1989) 1171; 

Z. Banach and S. Piekarsky, J. Math. Phys. 35 (1994) 

4809. 

S. Shahidi and H. R. Sepangi, Int. J. Mod. Phys. D 20 
(2011) 77, larXiv:1009.5458l 

T. Matos and F.S. Guzman, Class. Quant. Gravity, 18 

(2001) 5055, |arXiv:gr-qc/0108027| 

L. G. Cabral-Rosetti, T. Matos, D. Nunez and R. 
A. Sussman, Class. Quant. Gravity, 19 (2002) 3603, 
|arXiv:gr-qc/0112044| 

J.E. Lidsey, T. Matos and L. Arturo Urena-Lopez, Phys. 
Rev. D 66 (2002) 023514, [arXiv:astro-ph/0111292| 
T. Matos and D. Nunez , Rev. Mex. Fis. 51 (2005) 71, 
|arXiv:astro-ph/0303594| 



[32] K. Lake, Phys. Rev. Lett. 92 (2004) 051101, 
|arXiv:gr-qc/0302067| 

[33] A. D. D olgov and M. Kawasaki , Phys. Lett. B 573 
(2003) 1, |arXiv:astro-ph/0307285| 

[34] V. Far aoni and S. Nadeau Phys. Rev. D 72 (2005) 
124005, |arXiv:gr-qc/0511094| 

[35] T. Chiba, T. L. S mith and A. L. Erickcek Phys. Rev. D 
75 (2007) 124014, |arXiv:astro-ph/0611867| 



[30] 
[31] 



[36] V. Faraoni and N. Lanahan-Tremblay, Phys. Rev. D 77 
(2008) 108501. larX~iv:0712.3252l 

[37] S. Capozzie llo, V. F. Cardone and A . Troisi, JCAP 0608 
(2006) 001, |arXiv:astro-ph/0602349| 

[38] Y. Sobouti, Astronomy and Astrophysics 464 (2007) 



12 



921, |arXiv:astro-ph/060 3302 



[40] I. Shapiro, Phys. Rev. Lett. 13 (1964) 789. 



[39] S. Weinberg, Gravitation and Cosmology: Principles and 
Applications of the General Theory of Relativity, Wiley, 
New York, 1972. 



[41] J.Binney and S. Tremaine, Galactic Dynamics, Prince- 
ton University Press, 1987. 



